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Abstract 

Nuclear electromagnetic currents are derived in time-ordered perturbation theory within an 
effective-field-theory framework including explicit nucleons, A isobars, and pions up to one loop, or 
N 3 LO. The currents obtained at N 2 LO, i.e. ignoring loop corrections, are used in a study of neutron 
radiative captures on protons and deuterons at thermal energies, and of A=2 and 3 nuclei magnetic 
moments. The wave functions for A=2 are derived from solutions of the Schrodinger equation 
with the Argonne v%s (AV18) or CD-Bonn (CDB) potentials, while those for A=3 are obtained 
with the hyperspherical-harmonics-expansion method from a realistic Hamiltonian including, in 
addition to the AV18 or CDB two-nucleon, also a three-nucleon potential. With the strengths 
of the A-excitation currents occurring at N 2 LO determined to reproduce the n-p cross section 
and isovector combination of the trinucleon magnetic moments, we find that the cross section and 
photon circular polarization parameter, measured in n-d and ft-d processes, are underpredicted by 
theory, for example the cross section by (11-38)% as the cutoff is increased from 500 to 800 MeV. 
A complete analysis of the results, in particular their large cutoff dependence, is presented. 

PACS numbers: 12.39.Fe, 13.40.-f, 25.10.+S, 25.40.Lw 



I. INTRODUCTION 



The present work is the first stage of a research program aimed at studying electromag- 
netic observables of light nuclei (A < 8), and particularly radiative capture processes in the 
three- and four-nucleon systems, within a theoretical approach in which many-body elec- 
tromagnetic current operators derived in chiral effective field theory (xEFT) [l], 0, HJ are 
used in transition matrix elements between nuclear wave functions obtained from realistic 
Hamiltonians with two- and three-body potentials. This "hybrid" approach was adopted 
in a recent calculation of the astrophysical factor for the p-p and p- 3 He fusion reactions by 
weak capture at the keV energies relevant in the interior of the Sun 

Neutron and proton radiative captures on 2 H, 3 H and 3 He are particularly challenging 
from the standpoint of nuclear few-body theory. This can be appreciated by comparing 
the measured values for the cross sections of thermal neutron radiative capture on 1 H, 
2 H, 3 He. Their respective values in mb are: (332.6 ± 0.7) [5], (0.508 ± 0.015) @, and 
(0.055±0.003) 0- Thus, in going from A=2 to 4 the cross section has dropped by almost four 
orders of magnitude. These processes are induced by magnetic dipole transitions between 
the initial two-cluster state in relative S-wave and the final bound state. The 3 H and 4 He 
wave functions, respectively ^3 and ^4, are approximately eigenfunctions of the magnetic 
dipole operator /x, namely fi z ^3 — H P ^3 and [1^4 ~ 0, where /i p =2.793 n.m. is the proton 
magnetic moment — the experimental value of the 3 H magnetic moment is 2.979 n.m, while 
4 He has no magnetic moment. These relations would be exact, if the 3 H and 4 He wave 
functions were to consist of the symmetric S-wave term only. In fact, tensor components 
in the nuclear potentials generate significant D-state admixtures, that partially spoil this 
eigenstate property. To the extent that it is approximately satisfied, though, the matrix 
elements (^3 | \i z | \I/i +2 ) and (\J/ 4 | \i z \ ^1+3) vanish due to orthogonality between the 
initial and final states. This orthogonality argument fails in the case of the deuteron, since 
then fi z ^2 — (/^p — Hn) (p2(S) X0V0 j where Xm s an d Vm t are two-nucleon spin and isospin 
states, respectively. The magnetic dipole operator can therefore connect the large S-wave 
component 02(S) of the deuteron to a T=l 1 So n-p state — the orthogonality between the 
latter and the deuteron follows from the orthogonality between their respective spin-isospin 
states. 

As a result of this suppression, the n-d, p-d, n- 3 He, and p- 3 H radiative (as well as p- 3 He 
weak) captures are very sensitive to small components in the wave functions, particularly 
the D-state admixtures generated by tensor forces, and to many-body terms in the electro- 
magnetic (and weak) current operators. 

There have been in the past several calculations of these processes in the conventional 
framework — referred to as the standard nuclear physics approach (SNPA) in Ref. 0] — see [s| 
and references therein. Some of these studies, in particular the recent ones of Ref. 0] , have 
used accurate (essentially exact) bound and continuum wave functions corresponding to re- 
alistic Hamiltonians, which provide an excellent description of A=3 and 4 binding energies 
and radii, as well as of a variety of low-energy scattering observables (see [Io| and refer- 
ences therein). The electromagnetic current operator includes, in addition to the standard 
convection and spin-magnetization terms of individual protons and neutrons, also two- and 
three-body terms, constructed from, respectively, the two- and three-nucleon potentials so 
as to satisfy exactly current conservation (CC) with them. The method by which this is 
achieved has been improved over the years fljj , and its latest implementation is discussed at 
length in Ref. 0] — for an alternative formulation, though, see Ref. fl2j ]. It is not unique, since 



obviously the CC relation puts no constraints on the transverse component of the current. 
Nevertheless, it does generate two- and three-body terms, whose behavior, particularly at 
short range, is consistent with that of the corresponding potentials. This behavior in the lat- 
ter is ultimately "determined" by reproducing a set of experimental two- and three-nucleon 
scattering data and binding energies. 

These currents have been shown to provide a very satisfactory description of a variety of 
electronuclear properties, including, for example, n-p capture [9[] and deuteron photodisin- 
tegration at low energy 13j, and magnetic moments of A=3-7 nuclei [14(1. Yet, they lead 
to ~ 10% overestimate of the experimental cross section in n-d capture [9|. The discrep- 
ancy between theory and experiment increases to ~ 60% in the case of the n- 3 He capture 
cross section [l5|, although this earlier study, in contrast to that of Ref. [§], is not based 
on wave functions derived from the latest generation of potentials. The one-body (impulse- 
approximation or IA) term of the magnetic dipole operator gives, respectively, only 45% 
and 10% of the n-d and n- 3 He cross-section experimental values because of the suppression 
mentioned above. 

Electromagnetic currents up to one loop corrections have been derived in xEFT within 
the heavy-baryon formalism by Park et al. in Ref . |16lj . These currents were used in hybrid 
calculations of the n-p capture cross section 

0, 

Il7l|. spin observables in n-p capture (l8| , 
and, more recently, magnetic moments of the deuteron and trinucleons 17 . 



In the present work, we derive the electromagnetic currents by including, in addition to 
nucleon and pion, also A-isobar degrees of freedom. Thus we assume that the A-nucleon 
mass difference, just as the pion mass, is of the same order as the low momentum scale 
generically indicated by Q. Formal expressions up to one loop are obtained in time-ordered 
perturbation theory by employing non-relativistic Hamiltonians derived from the chiral La- 
grangian formulation of Refs. [ITS S]- The present study is similar to that of Ref. 16 



albeit it uses a different formalism. Various aspects of the calculations are discussed in con- 
siderable detail. However, a discussion of renormalization is not given here: it is deferred 
to a later publication [l9j]. It is nonetheless opportune to comment on it. There are two 
stages of regularization necessary in the one loop calculations: the first is the usual regu- 
larization of the one-loop corrections to the potential and to the currents, and the second 
is the regularization necessary for solving the Schrodinger equation and for the calculation 
of the current matrix elements. All this must be followed by corresponding renormalization 
procedures. In the calculations to follow at next-to- next-to-leading order (N 2 LO), the only 
loop corrections needed are those of the one-body current, which can be absorbed into the 
proton and neutron magnetic moments and electromagnetic radii. Therefore the required 
regularization is in the calculation of the matrix elements of the two-body components of 
the current. The latter is implemented as usual through a short-range cutoff parameter and, 
although a full fledged renormalization is not carried out at this point, we expect that a 
choice of the cutoff equal to that in the input potentials will give realistic values for the 
matrix elements. 

These N 2 LO currents are used to calculate the magnetic moments of A=2 and 3 nuclei, 
and thermal neutron radiative captures on protons and deuterons. Realistic two- and three- 
nucleon (for A=3) potentials are used to generate the bound and continuum wave functions. 
To have an estimate of the model dependence arising from short-range phenomena, the 
variation of the predictions is studied as function of the short-range cutoff parameter men- 
tioned above, which is used to regularize the two-body operators, as well as of the input 
potentials — either the Argonne v is (AV18) or CD-Bonn (CDB) in combination with 



respectively the Urbana IX |22| or Urbana IX* [23j — used to generate the wave functions 
(the AV18 and CDB have rather different short-range behaviors). 

We find that the N 2 LO calculations do not provide a satisfactory description of the 
experimental data, particularly for the suppressed process 2 H(n, 7) 3 H. This clearly points to 
the need of including loop corrections. However, it remains an interesting question whether 
these corrections will resolve the present discrepancies between theory and experiment. 

This paper is organized into eight sections and four appendices. In Sec. [TT1 we list, after 
defining our notation and conventions, the relevant strong- and electromagnetic-interaction 
Hamiltonians, obtained from chiral Lagrangians with nucleons, A isobars, and pions, while 
in Sec. IIHI we derive the nuclear electromagnetic current up to N 2 LO, i.e. ignoring loop 
corrections, in momentum space, and give the configuration-space representation of its op- 
erators in Sec. HVl Section |V] consists of two subsections: the first contains a derivation of 
one-loop two-body currents, while the second lists the two-body currents at N 3 LO, implied 
by four-nucleon contact Lagrangians involving two gradients. In Sec. I VI I we show that the 
currents up to N 3 LO are conserved when used in combination with xEFT potential including 
corrections up to one loop. In Sec. IVII[ we present and discuss results for the magnetic mo- 
ments of the deuteron and trinucleons, and for the radiative captures of thermal neutrons on 
protons and deuterons. Finally, in Sec. IVIIII we summarize our conclusions and outline the 
next stage in the research program initiated here. A number of details are relegated in the 
Appendices, including: expressions for the vertices associated with the interaction Hamilto- 
nians of interest (Appendix [A]); a collection of formulae relevant for the configuration-space 
representation of the N 2 LO operators (Appendix [B]) ; a listing of the analytical expressions 
for the one- loop currents involving A-isobars in the intermediate states (Appendix [U]) ; and, 
lastly, a listing of the four-nucleon contact Hamiltonians (Appendix ID]) . 



II. PRELIMINARIES 



Before listing the interaction Hamiltonians, it is useful to define our notation and con- 
ventions. In the Schrodinger picture adopted in the present study, the isospin triplet of pion 
fields 7r a (x) and their canonical conjugates n a (x), with a — x,y, z, are represented as 
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where the annihilation and creation operators c Pja and cp a satisfy standard commutation 
relations, and uj p = (p 2 + m 2 ) 1 / 2 . The nucleon and A-isobar fields, respectively iV(x) and 
A(x), with their corresponding canonical conjugates iN*(pc) and zA'(x), are given, in the 
non-relativistic limit of interest here, as 

J2b Ptar e^ X ar , (2.3) 
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where the 6's and gPs are annihilation operators for nucleons and A isobars, and x<tt and 
Xu A t a denote their respective spin-isospin states. Again, the 6's and b'% and similarly the 
ofs and c^'s, satisfy the standard anticommutation relations, appropriate for fermionic fields. 



Standard time-ordered perturbation theory is used to calculate the transition amplitude: 
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where | NN; 7) (| N'N')) represents, in a schematic notation, an initial (final) state con- 
taining two nucleons and a photon (two nucleons only) of energy Ei (Ef = Ei), Hq is the 
Hamiltonian describing free pions, nucleons and A-isobars, and Hi is the Hamiltonian de- 
scribing interactions among these particles as well as their couplings to the electromagnetic 
field. The evaluation of this amplitude is carried out in practice by inserting complete sets of 
eigenstates of H between successive terms of Hi in the expansion above. However, since we 
are only interested in retaining irreducible contributions, we omit from these intermediate 
states those that contain only two nucleons (this aspect of the present calculations is elab- 
orated on in Sees. IIII Al and IV Al) . Thus, the nuclear electromagnetic current j of interest 
here is obtained from 

(N'N' I T I NN; 7) Ududbic = — • j , (2.6) 

where q, ui q , and e q ^ denote respectively the photon momentum, energy, and polarization, 
and only first-order contributions in the electromagnetic interaction are retained in the 
evaluation of the transition amplitude. 



A. Pion, nucleon, and A-isobar interaction Hamiltonians 



The interaction Hamiltonians involving nucleon, A-isobar, and pion fields are derived 
from the effective chiral Lagrangian approach formulated in Refs. [3, 0]. Their explicit 
expressions are given by 
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where qa — 1-25 and F n ~ 186 MeV are the nucleon axial coupling constant and pion decay 
amplitude, respectively, Ha is the 7^ A coupling constant, and S a and T a are transition spin 
and isospin operators, converting a nucleon into a A isobar and satisfying 
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and similarly for Tj"T{, (note that in Ref. [2] the isospin transition operator is half that 
defined here). 

In addition to these, there is a set of four-fermion contact interactions described by 
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and S/2 is the spin-3/2 operator. As it will become clear below, terms involving more than 
two A-isobars are not needed in the present study. Finally, when discussing the renormaliza- 
tion of the two-body currents at tree level in a later work [19J , we shall also need to consider 
the following Hamiltonians involving three- and four-pion interactions, 

^-|dx7r 2 (x)iVt( x ) [o-.V7r a (x)] r a iV(x) , 
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obtained by including corrections up to 7r 2 (x)/.F 2 in the expansion of D 1 factors, where 
D = 1 + 7r 2 (x)/F 2 , entering the chiral Lagrangians [l|. 

B. Electromagnetic interactions 

The charged pion field is defined as 



^t( x ) = -j= M x ) T i7r B (x)] , 

and minimal substitution, 

Vvr T (x) -> [V^eA(x)]7r T (x) , 
in the pion-derivative couplings leads to the interaction Hamiltonians: 
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where e(> 0) is the electric charge, and the transverse vector field A(x) (in Coulomb gauge) 
is expanded as 
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The linear polarization (unit) vectors e p l , e p 2 form along with p a right-handed orthonormal 
system of axes, e p l x e p 2 = p. 

The interactions of individual nucleons and A-isobars with the electromagnetic field are 
described by the following Hamiltonians: 
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where k,s and Ky are the isoscalar and isovector combinations of the anomalous magnetic 
moments of the proton and neutron («g = —0.12 n.m. and Ky = 3.706 n.m.), and /x* is 
the iVA-transition magnetic moment (fi* ~ 3 n.m.). In the 7A^A term we only take into 
account the dominant magnetic dipole (Ml) coupling, ignoring the much smaller Coulomb 
(C2) and electric quadrupole (E2) couplings. The expressions in Eqs. (12.281) and (12.291) result 
from considering the non-relativistic limit of the effective Hamiltonians with non-minimal 
couplings 
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where i^n^x) and ^ (x) are the spinor and spinor- vector fields describing the nucleon and A 
isobar, and F uX (x) is the electromagnetic field tensor. The Bjorken and Drell conventions 24 
are used for relativistic four- vectors, 7-matrices, and Dirac spin- 1/2 spinors, except that the 
latter are taken to be normalized as u'(p, s)u(p, s) = 1. The Rarita-Schwinger spin-3/2 
spinors are defined as 
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where in the particle rest-frame the four- vector & is space-like, e M = (0, e A ), and A = ±1,0 
denote spherical components. 



C. Power counting 



We denote generically by Q a "small momentum", i.e. Q <C M, where M ~ 1 GeV is 
the typical hadronic mass scale, and consider as effective degrees of freedom the nucleon, 



A-isobar, and pion. Thus, we assume that the pion mass and the mass difference between 
the A-isobar and nucleon are both of order Q, m n ~ Q and ira — ~ Q. However, 
the photon energy u q is assumed to be suppressed by an additional factor Q/M relative to 
this small-momentum scale, i.e. u q ~ Q 2 /M. A generic coupling constant of dimensions 
(energy)" is assumed to scale with M as g = g M a with the expectation that g — 1. 

Contributions to the transition amplitude in Eq. (12.61) can be organized as an expansion 
in powers of Q/M The power counting implied by the interaction Hamiltonians in the 
previous two sections can be easily inferred by examining the structure of their associated 
vertices, listed in Appendix |A] This power counting is summarized in Table [H 
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TABLE I: Powers of Q, the small momentum scale, associated with the vertices from the strong- 
and electromagnetic- interaction Hamiltonians of Sees. Ill AMI Bl 

In the perturbative series, Eq. (12.51) . a generic irreducible contribution will be charac- 
terized by a certain number, say N, of vertices, each scaling as Q ai x Q^ 13 ^ 2 . . . , N), 
where ai is the power counting in Tableland /3, is the number of pions in and/or out of 
the vertex (this last factor is associated with the 1 / y/2 uj p included in the pion field), a cor- 
responding N—l number of energy denominators, and possibly L loops. Each of the energy 
denominators will involve pion energies and/or AiV mass differences, both of order Q, as 
well as kinetic energies of nucleons and/or A isobars, which, however, are suppressed by a 
further Q/M factor relative to Q. Loops on the other hand will contribute a factor Q 3 each, 
since they involve integrations over intermediate three momenta. Hence the power counting 
associated with such a contribution is 

irreducible contribution = (j[Q ai ~ M2 ^ X x Q 3L . (2.34) 

When one expands the nucleon propagators in powers of the kinetic energy around the static 
limit, terms with higher powers of Q appear as the kinetic energy is order Q 2 . 

The power counting of Eq. (|2.34[) can also be obtained by considering the Feynman 
diagram where the loop integrals are carried out in four dimensions. The power of Q of an 
irreducible diagram is then given by4L — 2n T — njq + Y^iLi a i-> where is the number of 
pion propagators, and the number of nucleon propagators. This equation together with 
the topological relation n n + = L + N — 1 leads to the power counting of Eq. (12.341) . 

Finally, the transition amplitude in Eq. (12. 6p can be represented diagrammatically as in 
Fig. [TJ The disconnected contributions in panels a) and b) will each contain a 5-function 
in the initial and final three momenta of one of the two particles, for example panel a) 
ex 6(p' 2 — P2), and will therefore be enhanced by a factor Q~ 3 relative to the connected (and 



irreducible) contributions in panel c). The power counting of diagrams a) and b) is then 
eQ x Q~ 3 = eQ~ 2 . In fact, these diagrams are the leading contributions to the nuclear 
electromagnetic current. This fact certainly fits in well with what is known empirically 
about, for example, magnetic moments of nuclei or radiative captures, such as the 1 H(n, 7) 2 H 
process considered later in this work. 
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FIG. 1: Schematic representation of the disconnected, a) and b), and connected, c), contributions 
to the NN^f — > NN amplitude. Solid and wavy lines denote nucleons and photons, respectively. 



III. CURRENTS UP TO N 2 LO 

In this section we derive the nuclear electromagnetic currents up to next-to-next-to- 
leading order (N 2 LO), that is e Q°. The relevant contributions are illustrated in Figs. [2] and [31 
The expressions below are given in momentum space, configuration-space representations 
are discussed in Sec. IIVI The momenta are defined as 
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where pj and p- are the initial and final momenta of nucleon i. The leading order (LO), 
eQ~ 2 , is given by the one-body current, panel a) in Fig. [2J 
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where q is the photon momentum, q = kj. 

The contributions to next-to-leading order (NLO), eQ -1 , are represented by diagrams 
b)-c) and d)-i) in Fig. [2J A straightforward evaluation of these diagrams in the static limit 
leads to the expressions: 
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where the momenta transferred to nucleons 1 and 2 add up to q, k x + k 2 = q. 

At next-to-next-to-leading order (N 2 LO), eQ°, there are two distinct contributions, illus- 
trated in Fig. [3l The first is due to (Q/M) 2 corrections to the one-body current in Eq. (13.21) . 
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FIG. 2: The complete set of time-ordered diagrams contributing at LO and NLO. Nucleons, pions, 
and photons are denoted by solid, dashed, and wavy lines, respectively. 

These are easily derived from a non-relativistic expansion of Eq. (12.311) : 
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At this point we should comment on the one- loop corrections to the one-body current, see 
Fig. HI They occur at NLO and N 2 LO, and are absorbed into the anomalous magnetic 
moments and electromagnetic radii of the proton and neutron. 

The second class of N 2 LO contributions, represented by diagrams b)-i), involve A-isobar 
excitation, and therefore vertices from H^na, H^na, and Hct,2- In the static limit, we find: 

4 r 2)Z k 2 - (n x t 2 ) z ct x x k 2 xq + 1^2, (3.6) 
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where A is the A-N mass difference, A=m& — m^, and use has been made of the identities 
in Eq. (I2.10p to eliminate spin- and isospin-transition operators in favor of Pauli spin and 
isospin matrices. 
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FIG. 3: The complete set of time-ordered diagrams contributing at N 2 LO. The square represents 
the (Q/M) 2 correction to the LO one-body current, while A isobars are denoted by thick solid 
lines, otherwise the notation is as in Fig. [2l 



We conclude this section by noting that the expressions in Eqs. (I3.3l) -( l3~4l) and (13.61) are 
the well known pion seagull and in-flight, and A-excitation currents commonly used in the 
literature (see, for example, the classic work of Ref. (2o|). 



A. Recoil corrections: cancellations at N 2 LO 



In the present formulation based on time-independent perturbation theory, there are in 
principle the additional N 2 LO contributions represented by diagrams a) and b) in Fig. [51 
However, these are exactly canceled by recoil corrections, also entering at N 2 LO, to the 
reducible diagrams c)-f). For example, the contribution of diagrams c)+d) is given by 



c)+d) 



j^yjyv-) K-jyjv(l) K-iVJv(2) 
Ei-Ep -E^-uj g + ir) 



Ei — E p — E 2 - uj q - uj k + ii] 



Ei-Ex-E 2 - uj q - u k + irj 



(3.8) 



where the V(iys denote the vertices from the interaction Hamiltonians relative to nucleon 
i, and the labeling of momenta is as illustrated in the figure. The initial and final energies 
Ei and Ef (Ei = Ef) are Ei—E\ + E 2 + ui q and Ef=E{ + E^ while E p is the energy of the 



intermediate nucleoli. These energies are all suppressed by Q/M relative to u k ~ Q, and 
therefore the denominators in square brackets can be expanded as 



2 Ei — E p — E 2 — u q 

u k u 2 k 



so that the contribution of diagrams c+d) now reads 

1 



c) + d) = V. 



-yNN\ 



Ei-E p -Ei—u q +i7) 



1 + 



2u k 



where ^(k) is the static one-pion-exchange potential (OPEP), defined as 



«*(k) 



VirNN (1) KrJVAr(2) 



g\ (Ti ■ k <t 2 ■ k 

>2 ^ + m 2 



Tl ■ T 2 



(3.9) 



(3.10) 



(3.11) 



The second term on the r.h.s. exactly cancels the contribution due to (the irreducible) 
diagram a) in Fig. [5j as can be easily surmised by noting that, in the static limit, the two 
energy denominators occurring in this diagram are each given by (— l/u k ), and therefore 



a) 



V yNN (l)v"(k) 
2u k 



(3.12) 



V 



a) 



b) 



FIG. 4: NLO and N 2 LO corrections to the one-body current. Notation is as in Fig. [2 



The discussion above becomes more delicate when the intermediate states describe fully 
interacting rather than free particles. Let 5v 77 be the recoil correction to the static OPEP 
v 7r , and let \(p) denote a bound or continuum state corresponding to Hq + v n with energy E. 
To first order in dv 7 * , the perturbed state \ is 

\^) =W) + -6v* | <p) =\ <p)+ 1 6<p) , (3.13) 

and the matrix element of the current operator j between initial and final states | ipi) and 
\ijjf) can be expressed as 

<V>/ 1 j I = (<Pf\i\ f*) + (<pi I J LO I Stpi) + (6<p f | j LO | 

= (ifi f | j | + (ifi f | j LO - 1 -5v- + h.c. | , (3.14) 
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FIG. 5: Time-ordered diagrams illustrating the cancellation of the irreducible contributions a) and 
b) by the recoil corrections to the LO diagrams c)-f). Notation is as in Fig. [2j 



where we have dropped terms of order (Sip) 2 as well as corrections beyond LO to the current 
in the matrix elements between | fij) and | ^v 9 /,*)- m the analysis following Eq. (13.8j) . in 
which the nucleonic intermediate states are free particles, the recoil correction can be written 

as 

5v" | free = (E t -H )^ + h.c. , (3.15) 

where uj is the pion energy. In momentum space the expression for 5v n \f me coincides with 
that implied by Eq. (|3.10p . namely 

v*(k) k-(K!-K 2 ) 

dV (k,Ki,K 2 ) |free= ~ , (3.16) 

2 u)k rriN 

with Kj defined as in Eq. (13.11) . 

If we assume that the nucleonic intermediate states describe fully interacting particles, 
i.e. they are eigenstates of H + v n , then it is plausible that the correction 5v n should be 
expressed as 

6v* = (Ei -H Q -v K )^- + h.c. , (3.17) 

from which it follows that 

;LO , -NLO ( vT ' -LO , t, „ I \ / , i u ;LO 



(^/ I J I V>i> = (<Pf I j j - [— j LU + h.c.j | + (if f | — j LU + h.c. I <*) . (3.18) 

The last two terms exactly cancel the two-body current contribution represented in Fig. [5] 
panels a)-b), namely the terms in brackets. Thus, if OPEP is taken in the static limit, 
Eq. ( 13. lip , as is the case for the calculations reported below, then the contributions of 
diagrams a) and b) should not be retained, since they are canceled by recoil corrections to 
OPEP. 



IV. CURRENTS IN CONFIGURATION SPACE 



The calculations of electromagnetic observables reported in Sec. IVHI are carried out in 
configuration space, and hence configuration-space representations of the current operators 
are needed. Those of the one-body operators in Eqs. (13. 2p and (13. 5p . generically denoted as 
jW, follow from 

j (1) (q)=/ I e ^i«+n)/2 ^-^-n) 5(kx-q) j«(ki,Ki) , (4.1) 

Jki JKi 

while those for the two-body current operators j^ 2 -* are derived from 

j (2) (q)=/ / e ^-r le ,k 2 .r 2 ^ (ki + k2 _ q) j(2)( kl)k2 ) , (4.2) 



2 



where the momenta kj and Kj are defined as in Eq. ( 13. II) . and we have introduced the 
notation 

<1P 5(...) = (27r) 3 5(...). (4.3) 



(2tt) 3 ' 

Note that Kj — > — iV^5(r^— r*), i.e. the configuration-space representation of the momentum 
operator. The LO current is then found to have the standard expression associated with the 
nucleon's convection and spin-magnetization currents, 

Ja°(q) = ^—(*n,i f-»Vi , e^l + in NA v x x q e^ 1 ) + 1^2, (4.4) 

where [. . . , . . .], denotes the anticommutator — this notation is also used in Eq. (14. lip below. 

At NLO and N 2 LO, however, the operators have 1/r 2 and 1/r 3 singularities (r is the 
interparticle separation), which need to be regularized in order to avoid divergencies in the 
matrix elements of these operators between nuclear wave functions. We adopt a simple 
regularization procedure [2j, i.e. a momentum-space cutoff. While its precise functional 
form is arbitrary, the choice made here of a Gaussian cutoff function, 

C A (p) = e~^ 2 , (4.5) 

with the parameter A < M ~ 1 GeV, is merely dictated by convenience, since it leads to 
analytical expressions for the Fourier transforms below. It is expected that this arbitrariness 
be of little relevance, since the dependence of theoretical predictions on variations in the 
cutoff is (or should be) largely removed by a renormalization of the theory free parameters, 
which are fixed by reproducing a given set of observables. We shall return to this issue later, 
in Secs.EnandlVnU 

The two-body currents at NLO are obtained as 

j^ L c°(q) = ^(rix r 2 ) z <r x (a 2 ■ V) / A (r) + 1 # 2 , (4.6) 

7T 

jd-P(q) = fa x t 2 ) z e^- R [a, ■ (V + zq/2)] [<r 2 ■ (V - iq/2)} Vg A (r, q) , (4.7) 



where r and R denote the relative and center-of-mass position vectors, respectively r = rj— r 2 
and R = (ri + r 2 )/2, and the functions f\(r) and g\(^, q) are defined as 



A(r) 



p 2 + m 2 



#A(r 



q) = J + \xe~ ix * rl2 J e~ ip r 



C A (p) 



[p 2 + L 2 (g; x)] 



2 ' 



(4.8) 
(4.9) 



with 



L(q; x) = yjml + (1 - x 2 ) q 2 /A . (4.10) 

Note that standard Feynman parametrization techniques have been employed to express g\ 
in the form given above. We defer to Appendix [B] for a listing of the formulae resulting from 
application of the gradients to f\(r) and ^(r, q). 

Finally, at N 2 LO the one-body term, Eq. (13.51) . reads 



jf L °(q) 



16 m% 



ejv,i 



2i(-V 2 + g 2 /4) (2Vi-(t 1 x q) 



- iVi-q (q + 2 CTl x Vi) , e iqri 



i e 



16 m 3 N 



(2 Vi-ffiX q) q 2 /2 - 2 (Vi X q) crj • Vi , e iqri 



- Vi • q (4<ri x Vi + q) 
+ 1^2, (4.11) 



while the two-body terms are given by 

4 r 2i , V - (n X r 2 ) 2 ^ x vl <r 2 • V/ A (r) + 1 # 2 , (4.12) 



•N 2 LO/ „\ _ _£Afl ^ i q . ri 

Jb -^ (q, " 8 9 miV AF 2 qX 



;N 2 LO f a _ • ^ -Pt iq-n 



4 r 2jZ cr 2 - (n x r 2 ) z cr x x cr 2 J /i A (r) + 1^2 
where h\(r) is simply the Fourier transform of the Gaussian cutoff function, 



hJr) 



-(Ar/2) 2 



(4.13) 



(4.14) 



In the limit A — > oo, h\(r) reduces to S(r). In this limit, as discussed later in Sec. IVHt the 
magnetic dipole operator derived from j^_^°(q) gives no contribution to nuclear electromag- 
netic observables — this follows from the antisymmetry of two-nucleon states. Smearing the 
5-function as in Eq. (14.131) is effectively including corrections of higher order than N 2 LO. 
We shall return to this issue in Sec. IVHi 



V. BEYOND N 2 LO: LOOP CORRECTIONS 

At N 3 LO (eQ), there are four classes of contributions: i) one- loop two-body currents, 
ii) currents from four-nucleon contact interactions involving two gradients of the nucleon 
fields, iii) one-loop renormalization corrections to tree-level two-body currents, and iv) 
(Q/M) 2 relativistic corrections to the NLO currents resulting from the non-relativistic re- 
duction of the vertices. We now turn to a derivation of the contributions in the first two 



classes. Those in the last two will be derived in a later work 19 



A. One-loop two-body currents 



In this section we consider one-loop two-body currents. Those involving pions and nucle- 
ons only in the intermediate states are illustrated by the diagrams in Fig. [61 where we show 
only one among all possible time orderings. Referring to this figure, we find 
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FIG. 6: Diagrams illustrating one-loop two-body currents. Only one among the possible time 
orderings is shown. Notation is as in Fig. [2j 



type a) 



-2i 



eg 2 A 

Ft 



T2,z (gi x q 2 ) + (n x t 2 ) z q 2 
UJ\ uj 2 {uji + 0J 2 ) 



+ 1#2 , 



(5.1) 



where the q^'s and u>i = (qf + m 2 ) 1 / 2 denote the momenta (with the flow as indicated in the 
figure) and energies of the exchanged pions, and the integration is on any one of the q^'s, 
the remaining q^'s with j ^ i being fixed by momentum-conserving 5-functions. Type b) 
diagrams give 



type b) = 2i 



c9a 



qi - qs 



UJi + UJ 2 + ^3 



F£ J UJ\ L0 2 UJ 3 {UJi + U)2){U)\ + UJ3) (uj 2 + L0 3 ) 



(ti x t 2 ) z qi ■ q 2 



r 2 , z <T\ ■ (qi x q2 y 



+ 1#2 . 



(5.2) 



Next, the contributions of type c)-e) diagrams are: 



type c) = -i 



( T i x r i)z 



qi - q3 + uj 2 + uj 3 ) - 3 uji u 3 

UJ1UJ3 (UJ\ + U)2){U\ + UJ 3 )(ijJ 2 + UJ3) 



(5.3) 



type d)=-2i^/ 



O^i cu| (^1 + ^2) 



vTi x r 2 ) 2 q 2 (qi ■ q 2 ) + 2 r 2 , 2 qi ■ q 2 (<Tl x q 2 ) 



2ri j2 q 2 cr 2 • (q x x q 2 ) 



+ 1 



(5.4) 



type e) = 2 z 



e 9A 

F 4 



qi - q3)/(^i, w 2 , W3) (Ti X t 2 ) 2 (qi ■ q 2 )(q 2 • q 3 ) 



+ 2 r 2)2 (q 2 ■ q 3 ) a x ■ (q 2 x q x ) + 2 r M (qi • q 2 ) cr 2 • (q 3 x q 2 ) 



(5.5) 



where the function f(ui,U2, 003) containing the pion energy factors from field normalizations 
and energy denominators for diagrams of type e) is defined as 



f(u U UJ 2 ,UJ 3 ) 



U)\ LU 2 LU 3 (lUi + UJ 2 )(uJi + UJ 3 )(U 2 + LU 3 ) 



OJl UJ 2 + UJ 2 UJ 3 + U)\ UJ 3 
0J\ UJ 2 uj 3 



{uji + u 2 ) (u 2 + u 3 ) {uj\ + uf) oo 2 u)\ + u 2 + oo 3 

1 o o I" ' 



Ujf Ul 2 U) 3 



LOi U 3 



(5.6) 



Lastly, diagrams of type f ) and h) vanish, since the integrand (in the static limit) is an odd 
function of the loop momentum qx, 



type f) and h) oc / x (spin— isospin structure) . 
However, the contributions of type g) and i) diagrams read: 



type g) 



2 1 



eg 2 A C T 



Tl X T 2 



q : - q 2 uj[ + Ui u 2 + LO2 



0J\ + 0J 2 



(5.7) 



(Ti ■ q 2 )(cr 2 • qi) , (5.8) 



type i) 



£9a 



Ti X T 2 



qi - q 2 u)\ + u)i uj 2 + uo\ 



Cso-i ■ (qi x q 2 ) 



2F% " J ujfuj 2 u)\ + uj 2 

+ i C T {(Ti ■ q 2 ) (a 2 ■ qi) + ((Ti ■ qi)(<x 2 ■ q 2 ) + i a 2 ■ (qi X q 2 ) 



1 # 2 . (5.9) 



A few comments are now in order. Firstly, the evaluation of the current operators resulting 
from the diagrams of Fig. [6] panels d)-e) and g) is carried out by including the recoil correc- 
tions of order e Q to the reducible diagrams shown in Fig. [7] panels a)-c). As an example, we 
consider the irreducible and reducible contributions represented in Fig. (H]d) and Fig. [7] a), 
respectively. We follow the procedure adopted in Sec. IIII Al and expand, in the reducible 
diagrams, the energies of the intermediate nucleonic states, which are suppressed by a factor 




FIG. 7: Diagrams illustrating the reducible one-loop two-body currents. Only one among the 
possible time orderings is shown. Notation is as in Fig. [2j 



Q/M with respect to the pionic energies Ui ~ Q. Up to order eQ, the current operator j red 
associated with the reducible box diagrams then reads 

jred= / = g ^ 7 ■ j NL °(qi) 

J Ei-Ei- E 2 + it) 



2 V nNN (2, q 2 ) K-atjv(2, qi) q?) V^atMI, qi) , (5.10) 



where v n (qz) and j NLO (qi) are the OPEP and pion-seagull current operators defined in 
Eqs. (13. lip and (13.31) . respectively. The V(i, q.,) denotes the vertex from the interaction 
Hamiltonian relative to nucleon % and a pion with momentum q,,, and E^ is the initial 
energy of the system, while E% and E 2 are the energies of the intermediate nucleons. The 
first term of Eq. (15.101) is then embedded in the iterated solution of the Lippmann-Schwinger 
equation, whereas the second term due to recoil corrections is retained and added to the 
irreducible contribution, j irr , which is given by 

Kriw(2, q 2 ) Kr7vv(2, qi) K-jvat(1; Qi) V^nn^, Qi) 



UJi uj 2 (^i + ^2) 
,2 1 , ,2 



+ [2 Ul+UJ . 2 + Ul ^ V„ NN (2, qi ) V nNN (2, qa) V nNN (l, q 2 ) 1^(1, qx) . (5.11) 

The first term above comes from the irreducible direct diagrams (in which, with reference 
to Fig. E]d), pion 1 is absorbed before pion 2), while the second term is from the crossed 
diagrams (in which pion 1 is absorbed after pion 2). Equation (15.111) can be further simplified 
expressing the product V 7vNN (2,qi) V n . NN (2,q 2 ) as 

V wNN (2, qi) V^ NN {2, qa) = [V nNN (2, qi), V^ NN (2, q 2 )]_ + V nNN (2, q 2 ) V nNN (2, q x ) , (5.12) 

to obtain 



,2 1 , ,2 



KrV7v(2, q 2 ) V n NN{2, qi) K-jvjv(1, q 2 ) V^-na^I, qi) 



/ 2— ^ — ^^[KvAf(2, qi), KAfv(2,q 2 )]_V r 7r Afv(l,q2)U 77rA r A r(l,qi) .(5.13) 

J Ui U>2{ UJ l + ^2) 



The complete current of type d) is then 



type d) 



1 7 : r - [VnNN{*, qi, 

WiW 2 Wl + 0J 2 ) 



V^jvjv(2,q 2 )]_K r jvjv(l,q2)^jvjv(l,qi) -h.c. , 

(5.14) 

where the h.c. term corresponds to including the diagrams in which the photon hooks up 
to the pion with momentum q 2 . Note that the recoil corrections exactly cancel the first 
term of Eq. (15.131) . leaving the term proportional to the energy factor associated with the 
crossed diagrams only. We find it interesting that these cancellations are also obtained for 
the current of type e), which reads 



type e) = / 4/(o> l3 u 2 , UJ3) [K-aw(2, q 3 ), V kNN (2, q 2 )]_ 
x V^ NN (1, q 2 )V KNN (l, qi)V^ 7r7r (qi, q 3 ) - h.c. 



(5.15) 



and it is therefore tempting to conjecture that they persist at higher orders. However, this 
statement has not been proven. 

Secondly, we observe that diagrams of the type shown in Fig. [8] are suppressed by an 
extra power of Q relative to those considered in this section, i.e. they are of order eQ 2 . For 
example, the diagrams of type a) give rise to the following current operator 



type a) in Fig. 8 



-fki(2K! 



i(Ti x ki) 



qi • q2 



+ 1 



(5.16) 



where the momentum Kj is as given in Eq. (13. ip . while those of type b) vanish, since they 
are proportional (5 az r ljb + S bz r ha - 2 S ab t 1jZ ) e abc r 2 , c = (e zbc + e bzc ) r 1>b r 2>c . 

Lastly, as a check, we have re-derived the nucleon-nucleon potential at the one-loop level 
(both with and without the inclusion of explicit A-isobar degrees of freedom), and have 
explicitly verified that it is in agreement with that obtained in Refs. j2| and 0]. In particular, 
we note that if recoil corrections to the reducible diagrams, for example box diagrams, are 
retained along with the contributions of irreducible diagrams, the resulting potential is in 
agreement with that derived in Ref. |3j with the method of unitary transformations (in this 
respect, see Sec. IVip . 
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FIG. 8: Diagrams illustrating N 4 LO contributions not included in the present work. Only one 
among the possible time orderings is shown. Notation is as in Fig. [2j 



We conclude this section by showing in Figs. M and [TU] the one-loop two-body currents 
involving one- and two-A intermediate states. A listing of the explicit expressions is given 
in Appendix O 
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FIG. 9: Diagrams illustrating one-loop two-body currents with a single A isobar in the intermediate 
states. Only one among the possible time orderings is shown. Thin, thick, dashed, and wavy lines 
denote nucleons, A isobars, pions, and photons, respectively. 

B. Currents from four-nucleon contact interactions 



In this section we report the N 3 LO contributions to the current operator from the contact 
electromagnetic interactions of Eqs. (IDlip - flDl9l) . We find: 



•N 3 LO 
JCT 7 



-eei 2 (2C' 1 -C 2 ) K 2 + 4 C'^ x + % C' A (<n + <r 2 ) x k 2 + %G' h a x x k x 

- % C 6 tr 2 x ^ + 2 (2 C 7 - C 10 ) (K 2 ■ <r 2 ) <ti + 2 (2 C£ - C' n ) (K 2 • m) <r 2 

- 2 C( 3 [(Ki ■ a" i) <r 2 + (K x ■ <r 2 ) o-x] + 2 (2 - C 12 ) K 2 (<n • <r 2 ) 

- 4C r i 4 K 1 (o-i-ff 2 )l +1#2. (5.17) 



Again, we note that the nucleon-nucleon potential generated by the contact interactions of 
Eqs. (ID2j) - flD10l) is in agreement with that obtained in Ref. (3|. 



f- ki p 2 H 


-k 2 














2 


1 




? 



Pl 



/) 



P2 




m) 





FIG. 10: Diagrams illustrating one-loop two-body currents with two A isobars in the intermediate 
states. Only one among the possible time orderings is shown. Notation is as in Fig. [9j 

VI. CURRENT CONSERVATION UP TO N 3 LO 



The nuclear electromagnetic current operator is related to the Hamiltonian through the 
continuity equation, which in the momentum space reads 



qj 



Pi 



+ 



+ V12 , P 



(6.1) 



where [. . . , . . .]_ denotes the commutator, q is momentum transfer, and p is the charge 
operator given to LO, in the notation of Eq. (I2.30p . by 



p = e (ei + e 2 ) . 



(6.2) 



It is well known (and easily verified) that the LO and NLO currents satisfy the continuity 
equation with, respectively, the kinetic energy terms and the LO (Q°) contribution to the 
potential, i.e. OPEP. The N 2 LO currents involving A excitation are purely transverse, and 
therefore do not enter the continuity equation, while those arising from relativistic correc- 
tions to the LO one-body term require the inclusion of these corrections also in the charge 
operator, in order for the continuity equation to be satisfied. We will not discuss them 
further here. 
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FIG. 11: Diagrams illustrating the reducible, panel a), and irreducible, panel b), two-body box 
potential. Only one among the possible time orderings is shown. Notation is as in Fig. [2J 



As an internal check, we have explicitly verified that the N 3 LO current operators obtained 
in the previous section satisfy the continuity equation with the N 2 LO (Q 2 ) contributions to 
the potential, induced by the interaction Hamiltonians given in Sec. Ill Al and Appendix ID 1[ 
For the purpose of illustration, we give more details on the calculation carried out for the 
currents of type d) and e) of Fig. El The potential generated by the box diagrams shown in 
Fig. [11] is given by 



^box = 



9 A f ^2+^3+ ^2 ^3 



2 n ■ T 2 (q 2 ■ q3) 2 + 3 cri • (q 2 x q 3 ) cr 2 • (q 2 x q 3 ) 



(6.3) 



F*J w| wf (u 2 + uj 3 ) 

where q 2 + q 3 = k 2 or — ki. We note that the recoil corrections to the reducible box diagrams 
have been included, consistently with our treatment in the previous section. 
Evaluation of the commutator of v^, ox with the charge operator gives 

2^1 



. ^box , P I 



-i e- 



Ti X T 2 



bj\ + U)\ + UJ 2 UJ 3 



(q 2 ■ q 3 ) 2 + 1^2 



(6.4) 



F£ " '*"J ^2^1(^2 + ^3) 

where k x + k 2 = q and q 2 + q 3 = k 2 . On the other hand, the l.h.s. of Eq. (16.11) for the 
currents of type d) and e) in Fig. [6] reads 

2 9a [ + ul + u 2 ujz 



qjd 



— 1 e- 



F$ J ul u\ (u 2 + u 3 ) 



[tx x t 2 ) 2 q 2 (q 2 ■ q 3 ) + 2 r 2 * (q 2 • q 3 ) (<Ti x q 2 ) 



+ 2 t 1z q 2 cr 2 • (q 3 x q 2 ) 



q+ 1 



(6.5) 



qjc 



1 e- 



(Tl X t 2 ) z (qi • q 2 )(q 2 • q 3 ) 



2 9 A f u 2 + u 3 + u 2 ^3 (jJ 2 + U 2 + U>1 U! 2 
J [ wfwf (W 2 + W3) ^1^1(^1 + LU 2 ) 

+2 T 2>z (q 2 • q 3 ) a x ■ (q 2 + 2 r M (q x ■ q 2 ) cr 2 ■ (q 3 x q 2 ) 



(6.6) 



where qi = q 2 + ki, and the factor q ■ (qi — q 3 ) f(u\, u 2 , u 3 ) = (u\ — u 3 )f(uJi, u 2 , u 3 ) has 
been written as in the square brackets of the last equation. Combining Eqs. (I6.5l) - fl6.6p . we 
obtain 

2 9a [u 2 2 + uj 2 + u 2 u 3 



q • (jd + jc 



-i e- 



(Ti x t 2 ) z (q 2 • q 3 )' 



F% J luIujI(uj 2 + uj 3 ) 
+2 T 2;Z (q 2 • q 3 ) <T! • (q 2 x q 3 ) + 2 r M (q 2 • q 3 ) cr 2 • (q 3 X q 2 ) 



+ 1 



(6.7) 



and the last two terms of the previous equation vanish. This is easily seen by changing 
q 2 — > k 2 /2 + q 2 (implying q 3 = k 2 /2 — q 2 ), and observing that the integrands are odd under 
q 2 — > — q 2 . Hence we are left with the first term which is equal to Eq. (16.41) . showing that 
the continuity equation is indeed satisfied. 



Pi + ki p 2 + k 2 




FIG. 12: Diagrams illustrating the reducible, panel a), and irreducible, panel b), two-body one-loop 
contact potential. Only one among the possible time orderings is shown. Notation is as in Fig. [2j 

Two closing remarks are in order. First, ignoring the recoil corrections in both the 
potential and currents leads to a violation of the continuity equation. This remains valid 
also for the current of type g) and the potential t> ct , generated by the diagrams shown in 
Fig. [12] and given explicitly by 

V« = \^T l .T 2 CT l .CT 2 j C ±. (6.8) 

Second, in hybrid calculations, such as those reported below, current conservation is 
not strictly satisfied. Assuming, however, that differences between the xEFT and realistic 
potentials occur at orders higher than N 2 LO, the N 3 LO currents derived here are then 
approximately conserved. 



VII. ELECTROMAGNETIC OBSERVABLES AT N 2 LO IN ,4=2 3 SYSTEMS 



We present results obtained for a number of low-energy electromagnetic observables in 
the A=2 and 3 nuclei using the current operators derived at N 2 LO. In the tables to follow, 
we denote respectively with LO, NLO, N 2 LO-RC, N 2 LO-A, and N 2 LO-A c the contributions 
calculated with the one-body current, or impulse-approximation (IA), of Eq. (I4.4p . the one- 
pion exchange two-body current at tree level of Eqs. (I4.6l) - (l4.7p . the relativistic correction 
to the LO current of Eq. (14.111) . the single A-isobar excitation current of Eq. (14.121) . and 
lastly the two-body current of Eq. (14.131) . due to the contact (NN)(AN) interaction. In the 
long-wavelength limit of interest in the present work, the LO and N 2 LO-RC currents are 
completely determined by the experimental values of the proton and neutron magnetic mo- 
ments, respectively +2.793 and -1.913 in units of nuclear magnetons (mm.). The NLO cur- 
rent involves the combination gA/F n , for which we adopt the value (m^ q a / F n ) 2 / (4tt)=0.075 
as inferred from the Nijmegen analysis of nucleon-nucleon elastic scattering data [26(. The 
NA transition magnetic moment is taken to be /i*=3 n.m. from an analysis of 7A data at 



resonance 



271 ]. The coupling constant Ha in the N 2 LO-A term is fixed by reproducing the 



width of the A resonance, Ha = 2.191 and qa = 1-267, while the (unknown) coupling 
constant Dt in the N 2 LO-A c term is expressed as 



Dn 



f 



, 9 Ah / 
F3 



(7.1) 



and the parameter /' is determined as described below. 

In order to have a realistic estimate of the model dependence of the results, we use 
cutoff values A in the range 500-800 MeV and wave functions corresponding to two dif- 
ferent nuclear Hamiltonians. The wave functions for A=2 are derived from solutions of 
the Schrodinger equation with the Argonne vi$ (AV18) [20] or CD-Bonn (CDB) [2l| two- 
nucleon potentials, while those for A=3 are obtained with the hyperspherical harmonics 
(HH) expansion method (see Ref. 10J and references therein) from a Hamiltonian including, 
in addition to the AV18 or CDB two-nucleon, also a three-nucleon potential, the Urbana-IX 
(UIX) model @- The AV18/UIX and CDB/UIX* combinations both reproduce the exper- 
imental 3 H binding energy. The former also reproduces most of the measured low-energy 
N-d scattering observables [l(|, with the notable exception of the vector analyzing power 
in N-d elastic scattering. Unfortunately, HH continuum wave functions for the CDB /UIX* 
combination are not yet available. The UIX* model 23| is a slightly modified version of 



the original UIX [22] (in the UIX*, the parameter Uq of the central repulsive term has been 
reduced by the factor 0.812). 

We consider the following two- and three-nucleon observables: the 1 H(n, 7) 2 H cross sec- 
tion at thermal energies, the deuteron magnetic moment, the isoscalar and isovector com- 
binations of the trinucleon magnetic moments, the cross section and photon circular polar- 
ization parameter R c measured in the radiative capture of (polarized, in the case of R c ) 
neutrons on deuterons at thermal energies. At N 2 LO there are no three-body currents. We 
also observe that at this order the only isoscalar terms are from the (one-body) LO and 
N 2 LO-RC operators, which are independent of the cutoff A. In Tables [Til and II III we list 





AV18 


CDB 


LO 


0.8469 


0.8521 


N 2 LO-RC 


-0.0082 


-0.0080 



TABLE II: Contributions in n.m. to the deuteron magnetic moment, obtained with the AV18 and 
CDB potential models. The experimental value is 0.8574 n.m.. 





AV18/UIX 


CDB/UIX* 


LO 


+0.4104 


+0.4183 


N 2 LO-RC 


-0.0045 


-0.0052 



TABLE III: Contributions in n.m. to the isoscalar combination of the trinucleon magnetic moments, 
obtained with the AVT8/UIX and CDB/UIX* Hamiltonian models. The experimental value is 
0.4257 n.m.. 

their contributions to the deuteron magnetic moment and isoscalar combination of the 3 He 
and 3 H magnetic moments. The N 2 LO-RC correction is (in magnitude) about 1% of the 
LO contribution but of opposite sign, so that its inclusion increases the difference between 



the measured and calculated values. As a result, the experimental deuteron and trinucleon 
isoscalar magnetic moments are underpredicted by theory at the (1.6-2.1)% and (3.0-4.7)% 
levels, respectively, depending on whether the CDB and CDB/UIX* or AV18 and AV18/UIX 
combinations are adopted in the A=2 and A=3 calculations. However, a recent calculation 
of these same observables (l7j | . based on variational Monte Carlo (VMC) wave functions cor- 
responding to the AV18/UIX Hamiltonian, finds the magnitude of the N 2 LO-RC correction 
somewhat smaller in A=2 (-0.0069 n.m.) and significantly larger in A=3 (-0.012 n.m.) than 
obtained here. The magnetic moment operator is derived via 



x j(q) 



q=0 



and from Eq. (14. lip the N 2 LO-RC term follows as 

[ 



,,N 2 LO 



A 



=1 L 



\p1 j 



+ e Ni Pi x (cr; x pi 



(7.2) 



(7.3) 



where p« = — zVj and are the linear momentum and angular momentum operators of 
particle i, and [. . . , . . .]+ denotes the anticommutator. The expression for this correction is 
different from that given in Ref . 17J , which is, in turn, different from that listed in an earlier 
work [18| by some of the same authors of Ref. [l7j . When compared to Eq. (17.31) above, for 
example, the term with the anticommutator is missing in Ref. 18], since the authors of that 
work normalize the spinors as uu = 1 rather than v)u = 1 as in the present work. These 
differences might partly explain the different contributions calculated here and in Ref. 17 
for the N 2 LO-RC correction. 

Next, we examine the 1 H(n, 7) 2 H radiative capture at thermal neutron energies. Various 
aspects of the calculations, which will not be discussed here, were reviewed most recently 
in Ref. 28| . The calculated values for the cross section are listed in Table IIVI As remarked 







AV18 






CDB 




A (MeV) 


500 


600 


800 


500 


600 


800 


LO 


304.6 


304.6 


304.6 


306.6 


306.6 


306.6 


NLO 


319.1 


319.6 


319.9 


321.3 


321.7 


321.9 


N 2 LO-RC 


317.4 


317.9 


318.2 


319.9 


320.3 


320.5 


N 2 LO-A 


321.9 


323.8 


326.3 


323.8 


325.3 


327.1 



TABLE IV: Cumulative contributions in mb to the 1 H(n,7) 2 H cross section at thermal neutron 
energy, obtained with the AV18 and CDB potential models and cutoff values in the range 500-800 
MeV. The LO and N 2 LO-RC contributions are cutoff independent, while the matrix element of 
the N 2 LO-A c operator vanishes. The experimental value is 332.6(0.7) mb from Ref. [B]. 

earlier in Sec. IIVI we note that the N 2 LO-A c current does not contribute, since the associated 
magnetic moment operator, 

2e/i*D T A 



,,N 2 LO 



9rriN A 



TLz tTi+T. 



{n X Tj) z (Ti X (Tj h A (r 



10/ 



(7.4) 



i<j=i 



is readily seen to vanish when acting on the 1 So n-p state If P[p PZ, and PL denote 
respectively the space, spin, and isospin exchange operators, where 

1 + <T j ■ (Tj 



P°. 



(7.5) 



and similarly for P[j, then P^ P i 



y y v 



-1 for a two-nucleon state. The identity 



<T, X (Tj 



(7.6) 



and the analogous one for Tj x Tj allow one to express the magnetic dipole operator in 
Eq. CO} as 



,,N 2 LO 



2e/i* D T 



A 



A (7* 



2 (r M CT, + r i)2 <Tj) - (r i)2 - r. 



(7.7) 



When acting on a two-nucleon state of even relative orbital angular momentum, the square 
bracket in the equation above reduces to (11 Z +Tj 2 ) (c^+cr,,), and therefore vanishes, since this 
state will have either S=0 and T—l or S=l and T=0. Indeed, in the limit h^r^) — > 5(rjj) 
the operator in Eq. (17.4ft gives no contribution. It is in this sense that one can interpret 
contributions at finite A as representing corrections beyond N 2 LO. 

The cutoff dependence and the different short-range behaviors of the AV18 and CDB wave 
functions lead to a cross section of (324.5 ±2.6) mb. Thus, at N 2 LO the experimental value, 
(332.6 ±0.7) mb |5j, is underpredicted by roughly 2.5%. The LO and NLO contributions cal- 
culated here are in agreement with those obtained for the AV18 in Refs. [13] — which uses the 
same form for the cutoff function — up to tiny differences presumably due to numerics. The 
N 2 LO-RC contribution, however, is found to be considerably larger (in magnitude) here than 
in Ref. [17J], alth oug h both studies predict the same sign for it, opposite to the LO contribu- 
tion. Reference [17j] suggests that corrections from current operators at N 3 LO might resolve 
the present discrepancy between experiment and theory (at N 2 LO), and possibly reduce the 
model dependence in the latter. This issue will be investigated in future work. Instead, 
rather than using the A width and 'j-N data to fix the values of the coupling constants 
Ha and /i* entering the N 2 LO-A current, Eq. (14.121) . we determine the combination //* Ua 
by reproducing the n-p radiative capture cross section. In fact, we make the replacement 
H* fiA — ► / /•** h-A, take fi* = 3 n.m. and Ha = 2.191 gA (as before) and fix / accordingly. The 
latter is listed in Table |V] for the various combinations of potentials and cutoffs considered 



A (MeV) 


AV18 


CDB 


500 


3.352 


3.234 


600 


2.471 


2.447 


800 


1.772 


1.814 



TABLE V: The parameter / required to reproduce, for a given value of the cutoff A, the experimen- 
tal 1 H(n,7) 2 H cross section at thermal neutron energy with the AV18 or CDB potential models. 
See text for notation. 



in this work. The resulting N 2 LO contribution becomes then as large as the NLO — a some- 
what unsettling feature of the present procedure. At NLO there is a significant cancellation 
between the contributions of the seagull, Eq. (14.61) . and pion- in-flight, Eq. (14. 7J , currents. 
This destructive interference persists also for the three-nucleon observables considered next. 
Lastly, it is worth noting that conventional calculations of the 1 H(n, 7) 2 H cross section based 
on two-body currents constructed to satisfy current conservation with the potential used to 
generate the wave functions, accurately reproduce the measured value 

Results for the isovector combination of the trinucleon magnetic moments are presented 
in Table IVII Note that the row labeled N 2 LO-A lists the contributions obtained with 



the strength of the N 2 LO-A current determined as in Table |Vj Consequently, they are 
significantly larger (in magnitude) and much less cutoff-dependent than those at NLO. The 
NLO contribution calculated in Ref. 17 with VMC wave functions and a cutoff of 600 MeV 



AV18/UIX CDB/UIX* 



A (MeV) 


500 


600 


800 


500 


600 


800 


LO 


-2.159 


-2.159 


-2.159 


-2.180 


-2.180 


-2.180 


NLO 


-0.156 


-0.197 


-0.238 


-0.113 


-0.156 


-0.200 


N 2 LO-RC 


+0.029 


+0.029 


+0.029 


+0.024 


+0.024 


+0.024 


N 2 LO-A 


-0.258 


-0.253 


-0.250 


-0.205 


-0.202 


-0.200 


Sum 


-2.544 


-2.580 


-2.618 


-2.474 


-2.514 


-2.556 



TABLE VI: Contributions in units of n.m. to the isovector combination of the trinucleon magnetic 
moments, obtained with the AV18/UIX and CDB/UIX* Hamiltonian models and cutoff values 
in the range 500-800 MeV. The LO and N 2 LO-RC contributions are cutoff independent. The 
experimental value is -2.553 n.m.. 

is -0.205 n.m., which is 4% larger than obtained here. This is most likely due to differences 
in the wave functions (we note, incidentally, that VMC wave functions are less accurate 
than HH ones). We determine the strength of the N 2 LO-A c current to reproduce, for a 
given cutoff A and Hamiltonian model, the experimental isovector magnetic moment. The 
resulting values for the parameter /', defined in Eq. ( 17. ip . are listed in Table IVHl The violent 
change of /' as the cutoff A is increased is due to the fact that the short-range behavior of 
the N 2 LO-A c current is governed by a Gaussian of half- width 2/ A. 



A (MeV) 


AV18/UIX 


CDB/UIX* 


500 


-3.036 


-38.57 


600 


+11.86 


-25.25 


800 


+51.13 


+3.485 



TABLE VII: The parameter f required to reproduce, for a given value of the cutoff A, the ex- 
perimental isovector combination of the trinucleon magnetic moments with the AV18/UIX or 
CDB/UIX* Hamiltonian models. 



However, with the values of the parameters / and /' fixed as discussed above, the current 
up to N 2 LO is now completely determined. We can therefore use it to make predictions for 
the cross section ot and photon circular polarization parameter R c measured in the reaction 
2 H(n, 7) 3 H. At thermal energies this process proceeds through S-wave capture predominantly 
via magnetic dipole transitions from the initial doublet J =1/2 and quartet J =3/2 n-d 
scattering states. In addition, there is a small contribution due to an electric quadrupole 
transition from the initial quartet state. We adopt here the notation and conventions of 
Ref. HH and define 



m 2 2 = Mi 



01/21/2 



M° 3/2 3/2 



644 = E, 



3/2 3/2 



(7.8) 



where M[ SJ and Ef SJ are the reduced matrix elements (RME's) of the magnetic and electric 
multipole operators of order £, normalized as in Eq. (6.3) of Ref. [29]. In terms of these 



RME's, the capture total cross section is given by 



2 a q 3 



, 2 \l m 22 + m 44 + e 44 , , 



(7.9) 



where a = e 2 /(4ir) is the fine structure constant, v Te \ is the d-n relative velocity, q is the 
energy of the emitted 7 ray, and is the nucleon mass. Similarly, the circular polar- 
ization P r resulting from S-wave capture of a neutron polarized along the direction P„ is 



proportional to the parameter R c [29], i.e. P\ 



q, where 



(7/2)|m 44 


2 + V8Re(m 22 ml 4 ) + 


[5/2)\e u 


2 + ^ 


v/24Im(m 2 2e^ 4 ) - Im(m 44 e^ 4 ) 




m 22 


2 + 


777,44 


2 + 


e 44 


2 



(7-10) 

The predicted RME's are listed in Table IVnTt the cross section and parameter R c in Table HXl 
Note that only results corresponding to the AV18/UIX Hamiltonian model are available at 
this time for the reason explained earlier. At LO the quartet m 44 is, in absolute value, 







m 22 






777.44 






e 44 




A (MeV) 


500 


600 


800 


500 


600 


800 


500 


600 


800 


LO 


-10.6 


-10.6 


-10.6 


13.5 


13.5 


13.5 


-0.14 


-0.14 


-0.14 


LO+NLO 


-13.2 


-12.5 


-11.5 


13.3 


13.4 


13.5 


+0.02 


+0.02 


+0.02 


LO+ • • • +N 2 LO-RC 


-12.0 


-11.3 


-10.3 


13.5 


13.6 


13.7 


+0.02 


+0.02 


+0.02 


LO+ • • • +N 2 LO-A 


-20.2 


-19.4 


-18.4 


12.6 


12.7 


12.8 


+0.02 


+0.02 


+0.02 


LO+ • • • +N 2 LO-A c 


-20.6 


-18.3 


-15.6 


12.6 


12.6 


12.7 


+0.02 


+0.02 


+0.02 



TABLE VIII: Cumulative contributions (in fm 3//2 ) to the reduced matrix elements (RME's) of the 
2 H(n, 7) 3 H reaction at thermal energies, obtained with the AVT8/UIX Hamiltonian model and 
cutoff values in the range 500-800 MeV. See text for notation. The 77122 and 777.44 RME's are purely 
imaginary, while the e 44 RME is purely real. 



about 27% larger than the doublet 77722- However, the contributions at NLO and N 2 LO are 
large and interfere constructively with those at LO for 777.22, while they are much smaller 
and interfere destructively for 77744. Consequently, the doublet m 2 2 at N 2 LO is found to be 
larger than the quartet m 44 by 63%-23% as the cutoff A is increased from 500 MeV to 800 
MeV. The calculation of the RME's is carried out with the Monte Carlo (MC) integration 



techniques of Ref. 29j, and the results reported in Table IVnTl are obtained from a random 
walk consisting of a large number (of the order of 2M) configurations. The statistical errors 
associated with these MC integrations are typically less than 2% for 77722 and much less than 
1% for 777 44 . However, they are 25% for e 44 at LO, and indeed much larger than the central 
value at N 2 LO, so that beyond LO the value of this RME is consistent with zero. We note 
that in S-wave capture the e 44 RME is predominantly due to transitions S( 2 H) — > D( 3 H) and 
D( 2 H) — > S( 3 H), where S and D denote S- and D-wave components in the 2 H and 3 H ground 
states. In the case of the AV18/UIX Hamiltonian, the contributions associated with these 
transitions interfere destructively, thus producing a small e 44 . This cancellation was found 
to be significantly model dependent in Ref. [29J, and use of CDB/UIX* wave functions would 
presumably produce somewhat different results for this RME in view of the considerably 
weaker tensor force of CDB relative to AV18 at intermediate and short range. 













R c 




A (MeV) 


500 


600 


800 


500 


600 


800 


LO 


0.229 


0.229 


0.229 


-0.060 


-0.060 


-0.060 


LO+NLO 


0.272 


0.260 


0.243 


-0.218 


-0.182 


-0.123 


LO+ • • • +N 2 LO-RC 


0.252 


0.241 


0.226 


-0.152 


-0.109 


-0.041 


LO+ • • • +N 2 LO-A 


0.438 


0.416 


0.389 


-0.432 


-0.418 


-0.397 


LO+ • • • +N 2 LO-A c 


0.450 


0.382 


0.315 


-0.437 


-0.398 


-0.331 



TABLE IX: Cumulative contributions to the cross section gt (in mb) and photon polarization 
parameter R c of the reaction 2 H(n,7) 3 H at thermal energies, obtained with the AYT8/UIX Hamil- 
tonian model and cutoff values in the range 500-800 MeV. The experimental values for ot and R c 
are respectively (0.508 ± 0.015) mb from Ref. 0] and -0.42 ± 0.03 from Ref. 0. 



At N 2 LO the cross section is underpredicted by theory by (11-38)% as the cutoff is 
increased from 500 MeV to 800 MeV. This rather drastic cutoff dependence is mostly due 
to the contribution of the N 2 LO-A c current. Indeed removing it leads to a much weaker 
variation of the cross section — roughly ±5% about the value obtained with A = 600 MeV 
(next to last row of Table IIXI) . It will be interesting to see to what extent, if any, loop 
corrections at N 3 LO will improve the present predictions, and in particular reduce the cutoff 
dependence. 

The photon polarization parameter is very sensitive to contributions of NLO and N 2 LO 
currents, which produce more than a sixfold increase, in absolute value, of the LO result, and 
bring it into much closer agreement with the measured value. All results listed in Table IIXI 
for R c (and <jt) include the small RME, although it only has a significant effect for 
the LO prediction (i? c =-0.060 versus -0.072 depending on whether is retained or not). 
The cause of the R c sensitivity to corrections beyond LO becomes clear by examining the 
expression given in Eq. (I7.10p . Neglecting e 44 , it reads 



1 + 3R C 



7/2 - \m 22 /m 4A \ 
1 + \m 2 2/m AA \ 2 



(7.11) 



and therefore the value attained by R c is driven by the ratio Itt^/w^I, which is about 0.79 
at LO and ranges from 1.63 to 1.23 for A = 500-800 MeV. 

We conclude this section by remarking that recent calculations of n-d capture observ- 



ables [3l|], based on an effective field theory formulated in terms of nucleon, deuteron, and 
triton fields with gradient couplings, seem to lead to predictions which are in agreement 
with data. However, it should be stressed that such a theory cannot be applied to other 
processes, for example the n- 3 He capture, without including additional degrees of freedom. 
It is in this sense of a more limited scope than the approach advocated in the present work. 



VIII. CONCLUSIONS AND OUTLOOK 

The goals of the present work were twofold: firstly, to derive the nuclear electromagnetic 
current up to one loop, i.e. up to N 3 LO, within an effective-field-theory approach including 
explicit nucleons, A-isobars, and pions, and secondly to implement this formalism in the 
calculation of a number of few-nucleon electromagnetic observables at low energy by using 



accurate wave functions — the so-called hybrid approach, advocated, for example, in Refs. [1, 
161 ]. This last objective has been partially carried out here, since results have been reported 
only at N 2 LO, i.e. ignoring loop corrections. 

Up to this order, the only isoscalar terms are those generated in a non-relativistic expan- 
sion of the one-body current: they provide a (cutoff-independent) 1% correction — relative to 
LO — to the deuteron and isoscalar combination of the trinucleon magnetic moments. This 
correction is of opposite sign to the LO contribution, and therefore increases the underpre- 
diction of the corresponding experimental values from (0.9 ± 0.3)% for the deuteron and 
(2.7 ± 0.9)% for the trinucleons at LO to, respectively, (1.9 ± 0.3)% and (3.8 ± 0.8)% at 
N 2 LO. The spread reflects differences in the short-range behavior of the AV18 and CDB 
potentials, in particular the weaker tensor components of the latter relative to the former in 
this range. 

At NLO, isovector terms arise from the pion seagull and in-flight contributions, while at 
N 2 LO, in addition to the relativistic corrections mentioned above, isovector terms due to 
A-isobar excitation are also obtained. The value for the combination of coupling constants 
{9a/Fit) 2 entering the NLO two-body currents is that inferred by an analysis of nucleon- 
nucleon elastic scattering data [26j. However, the strengths of the N 2 LO two-body A- 
excitation currents, i.e. the combinations /j,* Qa^aI an d /■** D T in Eqs. (14.121) and (14.131) . 
respectively, have been determined, as functions of the cutoff A and for the Hamiltonian 
model of interest, by reproducing the cross section for n-p capture and the isovector com- 
bination of the trinucleon magnetic moments. This current has then been used to make 
predictions — with the AV18/UIX model only — for the cross section ot and photon circular 
polarization parameter R c . The experimental a? ( | -R c I ) is found to be underestimated by 
11% (overestimated by 4%) for A=500 MeV and 38% (underestimated by 21%) for A=800 
MeV. We note that the parameter R c is mostly sensitive to the ratio of doublet to quartet 
magnetic dipole transition matrix elements 1 77122/^44] (the cross section is proportional to 
|m 2 2| 2 + \m AA \ 2 ). 

The results display a significant cutoff dependence, particularly so for the N 2 LO contri- 
butions associated with A isobar degrees of freedom. Indeed these contributions are much 
larger than those at NLO. This is partly due to the fact that the two NLO (pion seagull 
and in-flight) terms interfere destructively. For example, the seagull (in-flight) contributions 
to m.22 and 77744, i n units of fm 3 / 2 and for A=500 MeV, are respectively — 9.1i (+6.5 i) and 
— 0.8 i (+0.6 i), which add up to the NLO values — 2.6 i and —0.2i from Table IVTHi As a 
result Ot = 0.425 mb and R c = —0.425 at LO+NLO (seagull only), which should be com- 
pared to oy = 0.272 mb and R c = —0.218 at LO+NLO (seagull+in-flight) from the second 
row of Table IIXI 

The relatively large A-excitation contributions also point to the need for including loop 
corrections at N 3 LO, which these N 2 LO currents, because of the procedure adopted here to 
determine their strength, are implicitly making up for. This is also evident by examining the 
results for the n-p capture cross section in Table HVl Had we chosen to fix the N-A transition 
axial coupling constant and magnetic moment fi* via, respectively, the A width and 7- A" 
resonance, the contribution of the N 2 LO-A current would have been considerably 
smaller than that at NLO, and would have fallen more in line with naive expectations. 

The N 3 LO corrections will presumably reduce the cutoff dependence in the n-d predic- 
tions, and hopefully bring theory into more satisfying agreement with experiment. For 
the time being, we only observe that calculations [9j based on the AV18/UIX Hamiltonian 
with leading two- and three-body currents constructed consistently — in the sense of being 



exactly conserved — with, respectively, the AV18 two-nucleon and UIX three-nucleon po- 
tentials overestimated both <jt and |i? c | by about 10% in the 2 H(n, 7) 3 H process at thermal 
energies, while at the same time providing an excellent description of cross section data for 
the p-d radiative capture in the energy range from a few to 80 keV, and, in particular, of 
the astrophysical factor at zero energy extrapolated from these data. 

Thus, as already emphasized in the Introduction, very low-energy radiative (and weak) 
capture reactions involving three- and four-body nuclei constitute a crucial testing ground 
for the models used to describe the ground- and scattering-state wave functions — and indi- 
rectly, the underlying interactions which generate these wave functions — and the many-body 
electroweak current operators [§]. 

The next stage in the research program we have undertaken is to incorporate the N 3 LO 
operators derived here into the calculations of the captures and magnetic moments involving 
light nuclei (with mass number A < 8), and indeed to extend these calculations to also 
include p-d capture at energies up to a few MeV's, and possibly four-nucleon processes, in 
particular 3 He(n, 7) 4 He at thermal energies. Of course, at N 3 LO three-body currents also 
occur, and will need to be derived. Work along these lines is being pursued vigorously. 
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APPENDIX A: VERTICES 

The interaction Hamiltonians in Sees. Ill Al and III Bl are assumed to be normal-ordered. 
Explicit expressions for the associated vertices are easily derived (these expressions include 
the 1 / y / 2uJfc~ factors from pion fields) : 

i) Strong-interaction vertices: 



(pa>xa;M I Httna I p,x) 





.9a cr ■ k 
Ma S ■ k 



(A2) 



(Al) 




(A3) 



(pi,xi;p 2 >x 2 1 h ct,i I Pi,xi;P2,X2> = 

(p',x';pa,Xa I h ct,2 I Pi,xi;P2,X2) = 
(p, x; pk> x'a I h ct,s I p, x; pa, xa> = 



C a Ti a • r 2a , 

a=S,T 

D T Ti • T 2 <Ti • S2 , 
X! C' a Fia ■ T' 2a , 

a=S,T 



(p'l,A> Xi,a5 P2,A> X 2 ,A I #CT,4 I Pi, Xl5 P2, X2) = £>T T l ' T 2 Si • S 2 , 

(pa, x a ; p', x' I ^ct, 5 I p, x; pa, xa) = £>t t 1 ■ t 2 s x • s 2 , 



(p', x'; ki, a; k 2 , 6; k 3 , c \ H 3n \ p, x) 



(cr • ki r a 5 bc 



F% y/8 uj kl uj k2 UJ k . 
+ cr • k 2 r b 5 ca + cr • k 3 r c 5 afe ) , 
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(ki, a; k 2 , b; k 3 , c; k 4 , d | # 4?r | 0) 



hb hd (hv k 2 + hfiK + m l) 



Fl y/V6 uj kl u k2 u k3 u k4 
+ S ac hd (hfikz + k 2ll k% + ml 

+ Sad he (h^K + ^2^3 + m l) 

ii) Electromagnetic-interaction vertices: 



/ / / , 1 Tj 1 w 9 A cr e qA 

(p , x ; k, a I H ywNN I p, x; q, A) = e— — j== 7 =^ e zab r b , 



F n \^2u k v2w 5 
(P , X A ; k, a \ # 7 ^a I P, X5 q, A) = e— — j= j=- e zab T b , 



(k 1 ,a;k 2 ,b \ H 7wn | q, A) =ie 
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(p', X'; ki, a; k 2 , b \ H 1?TnNN | p, x\ q, A) 
e q A (p' + p) + i & x (p' - p) 



Fl ^Auj kl u k2 



2m N 



(hz n + 5 bz r a - 2 5 aft r 2 ) , (A14) 
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(ki, a; k 2 , b; k 3 , c; k 4 , d \ H^ 4n | q, A) = -ie 
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In these expressions p and pa denote nucleon and A-isobar momenta in spin-isospin states 
specified by \ and xa respectively, while the k's and a,b, . . . denote pion momenta in isospin 
states a,b, . . ., and q and A the photon momentum and polarization state. For brevity, on 
the r.h.s. of the equations above the spin-isospin states of the nucleon and A isobar as well 
as the ^-functions enforcing three-momentum conservation, are not shown explicitly In 
Eq. (jAlOp . the notation /cf/c JAt denotes the combination io^lO}.. — kj • k,. Finally, vertices 
involving A-isobar deexcitation into a nucleon are obtained by replacing S and T by their 
adjoint operators and T^, while vertices in which one or more pions are in the initial 
state are obtained from those listed in Eqs. (jXIJl -dAgj) . (jA9l> . fTAToD . (jAl3|l . and flAT6l) by 
replacing kj — > — kj and/or uo ki — > —cj^ (of course, the energy replacements are not to be 
carried out in the pion-field normalization factors). For example, 



(p',x';ki,a | H, 



nnNN 
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F 2 V 4 Vki 



■ €abcT~c 



(A19) 



APPENDIX B: CONFIGURATION-SPACE REPRESENTATION 



We list here the configuration-space representation of the two-body currents at NLO and 
N 2 LO. To this end, it is convenient to define z n = m^r, z A = Ar, zl = rL(q;x), 



z± = m w /A ± z A /2 , z± = L(q; x)/A ± z A /2 
and the complement error function 

^(z) = -^=l die"* 



(Bl) 



(B2) 



The dependence of zl and z± upon the variable x is understood. Then the functions f\(r) 
Eq. (14.81) . along with its first derivative f A (r), and Vg A (r,q), Eq. (14.91) . are given by 
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where only the transverse part of Vg A (r,q) (orthogonal to the photon momentum q) is of 
interest, and 



E q (x,r) 



8nz 2 A 
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In the limit A — > oo these functions reduce to: 

f (r) = ?±— f (r) = - r ^ — (l + z) 



. 7T 



(B7) 



The complete NLO current — the sum of the two contributions in Eqs. (14.61) and (14.71) — is 
then written as 



j NLO (q) , = e^(n x r 2 U e^f A (r) <t x {ct 2 ■ f) + e***f' A (r) o 2 (oi ■ f ) 
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where r = r x — r 2 and R = (1*1 + r 2 )/2, and the functions ^ with z = 1, . . . , 5 are defined 

as 

q) = J + l da;e-"' q r/2 ^1 - r-^j E, 
:,q) = ^J^dxe-^ 2 (1 + x) E q (x, 

\f \xe-^ (1 + x) {l-r^j E, 
:,q) = Iy +1 dse- i:rqr/2 (l - x 2 ) E g (x,r) 
#(r,q) = J^dxe-^ (3-3rA + r 2 -^) £, 
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The configuration-space representation of the N 2 LO current in Eq. (14.121) reads 
. e/i* g A h A 
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where /ig(r) = f\(r)/r and 
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APPENDIX C: ONE-LOOP TWO-BODY CURRENTS WITH A ISOBARS 



We begin by including a single A isobar in the intermediate states. The relevant diagrams 
are shown in Fig. [91 We find for type a) and b) diagrams: 



h 2 
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(Si • q 2 ) Si 
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V 
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where S and T are the spin- and isospin-transition operators defined in Eq. (12.101) . and A 
denotes — m^. The spin-isospin structures can be further simplified and expressed in 
terms of the Pauli matrices er and r. 

The contributions of type c)-e) diagrams can be written as 
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type d) = 
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where we have defined the one-pion-exchange transition potential NN — > AiV as 

^ (Si ■ q 2 ) (cr 2 • q 2 ) 



^ATViq2j 



Ti-r 2 
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the transition currents ^yNN — > AiV associated with the seagull 77riVA-coupling and pion- 
in-flight terms respectively as 



JA7v(qiy 



&(qa>qs) 



2 e- 



• 9a h A Si o- 2 ■ qi N 

F2 g 2 + m 2 

9a hA qi - q3 



(C7) 



(Si • qi) (a 2 • q 3 ) (T x x t 2 ) z , (C* 



and vJf A and j^r A are obtained from v\ N and j AAr by the replacements S # cr and T 
The contributions labeled jjf^ e are given by 
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where the functions f± (uix, uj 2 ) and / 2 A (co>i, u 2 , cu 3 ) denote the following combinations of pion 
energies and AN mass differences: 
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FIG. 13: Direct box diagrams with a single A isobar in the intermediate states. See text for 
discussion. Dashed lines with a cross represent v an transition potentials, or j AN and j A w N transition 
currents. 

At this stage, it is useful to comment on the structure of the contributions in Eqs. (1C3|) - 
( 1C5|) . The first terms on the r.h.s. of each of these equations are represented by the diagrams 
in Fig. [T3J i.e. (from top to bottom) an instantaneous interaction mediated by a transition 
potential, a A-nucleon energy denominator taken in the static limit, and a two-body current 
inducing excitation of a single A. They have a simple interpretation [15(|: for example, the 
type c) in Fig. O is the matrix element (AN \ j n &N I NN) evaluated by treating the final 
I AN) state in first order perturbation theory, 



(A' N' I v AN 
AN) =\ NN) + £ I A'N') { 1 



NN) 
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J AN 
A 
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One additional feature of these terms is that their configuration-space representations 
are particularly simple, since they are given, for example for type c), by the product 



—v\ N (r) j^ JV (q)/A, and similarly for type d) and e). Here v\ N (r) and j*Ajv(q) are the 
configuration-space representations of, respectively, the transition potential and the current 
in Eq. flCTD . 

The remaining terms j£ j e represent corrections to this picture, arising from crossed-box 
diagrams. In particular, the functions (o^o^) and f 2 (uj±, ^2, ^3) give, up to pion-energy 
factors from field normalizations, the sum of the energy denominators for the six and thirty 
crossed-box diagrams from, respectively, type c)-d) and e) contributions, 
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1 



2 uj\ to 2 
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[sum of 6 type c) or d) crossed — box diagrams] , (C15) 
[sum of 30 type e) crossed— box diagrams] . (C16) 
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Lastly, the contributions of type f )-i) diagrams vanish, while those of type j)-k) are written 
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where the contact (momentum-independent) transition potential v\ N has been defined as 

v c AN = D T T 1 - r 2 Si ■ a 2 , (C19) 

and the current jj - ' 1 is given by 
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The next set of contributions we consider includes two A isobars in the intermediate states, 
the relevant diagrams are displayed in Fig. [TUJ In analogy to Eqs. (1C3|) - (1C5I) . we write: 
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where the transition potential v AA and currents j AA and are obtained from Eqs. (1C6j) - 
(1C8j) by the replacements qa — > ^a and cr 2 (t 2 ) — > S 2 (T 2 ), and 
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sum of 60 type n) direct and crossed— box diagrams] ,(C26) 



while / A (o;i,u; 2 ) and / A (uii, lu 2 , ^3) as in Eqs. f)C15j) and (]C16j) . but for diagrams of type 1) 
and m), respectively. They are explicitly given by 
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We observe that, in contrast to the case of a single A, the energy denominators for the direct 
and crossed box diagrams do not add up to 1/(2 A) — as one would have naively expected— 
indeed, that is the reason for including this term within the square brackets of Eqs. (1C23j) 
and (1U24D . 

Next, we consider diagrams n)-s) in Fig. [TUJ which involve contact terms with two A's. 
The contributions of type n)-p) vanish, while those of type q)-s) read: 
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APPENDIX D: CURRENTS FROM CONTACT INTERACTIONS 



In this Appendix we list the four-nucleon contact interaction Hamiltonians involving two 
gradients of the nucleon fields. Minimal substitution, 

ViV(x) -> [V - i e eA rA(x)] N(x) , (Dl) 

in the nucleon-derivative couplings then leads to the corresponding electromagnetic- 
interaction Hamiltonians, which are listed as well. In the last section of this Appendix 
we report the expressions for the vertices induced by these Hamiltonians. 



1. Four-nucleon contact interaction Hamiltonians 



The four-nucleon contact interaction Hamiltonians with two gradients acting on the nu- 
cleon fields have the following expressions 0, 0] : 
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iV t (x)VA^(x)l • [[VAT(x)] f x <riV(x) 



(D2) 
(D3) 
(D4) 



i C 4 I dx 



-^CT2D,5 

-f^CT2D,6 
-^CT2D,7 



+ [[ViV(x)] f iV(x)] • [iVt( x )o- x ViV(x) 

iC' 5 Jdx [iV t (x)A^(x)] [[VA^x)] 1 • a x ViV(x) 

iC' 6 fdx [iV t (x)criV(x)] • [[VA^x)] 1 x ViV(x)] , 
(C' 7 Sik5ji + C'gSuSkj + CgSijSki) 

J dx [iV t (x) ( r fe 9 i A^(x)] [iV t (x)(j^ i iV(x)" 

+ [[9,iV(x)] t ( T fc iV(x)l [[9 j 7V( x )] t ( 7 i 7V(x) 



(D5) 
(D6) 
(D7) 



-^CT2D,8 — (C[ Sik5ji + C^^j^fej + C' 12 S i jS i 



Id J 



J dx [iV t (x) < r fc 9 i iV(x)] [[<9,iV(x)]UiV(x 



(D8) 
(D9) 



H, 



CT2D,9 



2 {dik^ji + 5ii$kj) + C[ 4 5ij8ki 



|dx [[^iV(x)] t ( 7 fe ^iV(x) + [^iV(x)] t ( 7 fc a t iV(x)] [jV t (x) < r,JV(x)] .(D10) 



2. Contact electromagnetic-interaction Hamiltonians 



Minimal substitution leads to the following contact electromagnetic-interaction Hamilto- 
nians: 



H, 



CT 7 ,1 



= — i e 



#CT 7 ,3 
-^CT 7 ,4 



H, 



CT 7 ,5 



Hi 
H, 



CT 7 ,6 - 
CT 7 ,7 = 



H, 



CT 7 ,8 



y dxA(x) • [iVt(x) (\f - V) iV(x)] [iV t (x)e J vAr(x)' 

+ [Art( x ) ejv iV(x)] [iVt(x) - V) iV(x) 

-ieC' 2 J dxA(x) • [iV t (x)e 7V iV(x)] [[ViV(x)] t iV(x) 

- [A^ t (x)ViV(x)] [A^ t (x)e JV Ar(x 

ieC^/ dxA(x) • [2A^ t (x)iV(x)] [jV f (x) (\f - V) e w iV(x) 

eC^ dxA(x) • [iVt(x) (\f + V) 7V(x)] x [n*(x)<t e N N(x) 

+ [iV f (x) e N N(x)] [iV f (x) (\f + V) x crjV(x) 

eC^ dxA(x) • [iV t (x)iV(x)] [jV f (x) (\f + V) x o-e w iV(x) 

eC^ dxA(x) • [iV t (x)o-iV(x)] x [iV f (x) (\f + V) e w iV(x)' 
-ie (C' 7 5 ik 5ji + C' s 5 u 5 kj + Cg5ij5 k i) 

J dx [^(x) [iVt(x) (d* - ^) a fe 7V(x)] [iVt( x )^ ejv iV(x) 

+A,(x) [iVt(x)a fcejv iV(x)] [iVt( x ) (fif - ^) ^TV(x) 



(Dll) 

(D12) 
(D13) 

, (D14) 
, (D15) 
, (D16) 



, (D17) 



dx 



^(x) 



-A: fx) 



iVt(x) < 7 fc 9 i JV(x)] [iVt(x)a ; e w iV(x) 
'iVt(x)t7 fcejv JV(x)l [[^■iV(x)] t ^iV( 



(D18) 



CT 7 ,9 



J dx L(x) [iVt(x) (fif - <%) ^ ejv iV(x)] [jV t (x)a l JV(x) 
+A(x) [iVt(x) (flj - ^) (7 fcej viV(x)] [iVt(x)^iV(x; 



(D19) 



3. Contact interaction vertices 



The vertices induced by the contact electromagnetic-interaction Hamiltonians are listed 
below. The notation is the same as in Appendix |Aj but for 



and 



(H C T 7 ,i) = (pi, xi; P2, x' 2 1 H C T~,,i I pi, xi; P2, X2; q, A) 



^ = 1, 



9 



(ffcr 7 ,i> = 2eC( 



(if, 



CT 7 ,2 



) = -e^ 



(#CT 7 ,3) 



2eC 3 



ei (p 2 + P 2 ) + e 2 (Pi + Pi) 
ei (P2 + P 2 ) + e 2 (Pi + Pi) 
ei (Pi + Pi) + e 2 (p 2 + Pa) 



e q A 

y2uL 



(Hct-,,4) = -ieC'4 (o"i + 0*2) x 



(#CT 7 ,5) 



-ieC' 



ieC' 



ei (P2 - p' 2 ) + e 2 (Pi - Pi 
ei <Ti x (pj - pi) + e 2 cr 2 x (p 2 - p 2 ) 
ei cr 2 x (pj - pi) + e 2 o"! x (p 2 - p' 2 ) 



e q A 

»/2u;. 



e q A 
e q A 



<#ct 7 ,7> = 2e 



C 7 [ei (P2 + P2) ■ cr 2 cri + e 2 (pi + pi) • <y 1 <x 2 ] 
Cs[ e i (P2 + P2) • CT i CT 2 + e 2 (pi + pi) • cr 2 cti] 



+C9 a x ■ a 2 [e 1 (p 2 + p' 2 ) + e 2 (pi + pi)] 



e q A 



c 'wl e i (P2 + P 2 ) • °"2 cr x + e 2 (pi + pi) • CT\ <T 2 ] 
-C^fei (p 2 + p' 2 ) ■ cri cr 2 + e 2 (pi + pi) ■ cr 2 o"i] 
-C( 2 <Ti ■ cr 2 [ei (p 2 + p' 2 ) + e 2 (pi + pi)] 



e q A 



C' 13 [e! (pi + pi) • cri cr 2 + ei (pi + pi) • cr 2 <xi 
+e 2 (p 2 + p' 2 ) ■ (Ti cr 2 + e 2 (p 2 + p' 2 ) ■ cr 2 trj 



(D20) 

(D21) 
(D22) 
(D23) 
(D24) 
(D25) 
(D26) 



(D27) 



(D28) 



+2 C' u <X! ■ cr 2 [ei (pi + pi) + e 2 (p 2 + P 2 ) 



e q A 



(D29) 
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